Magnetic field switching in parallel quantum dots 
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We show that the Coulomb blockade in parallel dots pierced by magnetic flux $ completely blocks 
the resonant current for any value of <£> except for integer multiples of the flux quantum $o- This 
non-analytic (switching) dependence of the current on "!> arises only when the dot states that carry 
the current are of the same energy. The time needed to reach the steady state, however, diverges 
when $ — > n$o. 

PACS numbers: 73.23.-b, 73.23.Hk, 05.60.Gg 



The system of two quantum dots coupled in parallel 
to two reservoirs has attracted a great deal of attention 
as a realization of a mesoscopic Aharonov-Bohm inter- 
ferometer 0, 0, 0| . Indeed such a system pierced by an 
external magnetic field (Fig. Q]) is an interference device 
whose transmission can be tuned by varying the magnetic 
field. In the absence of the interdot electron-electron in- 
teraction, the interference effects in the resonant current 
through this system are quite transparent. This is not 
the case, however, for interacting electrons 0,0- 

Consider for instance a strong interdot electron- 
electron repulsion — a Coulomb blockade. While the two 
dots may be occupied simultaneously in the noninteract- 
ing model, the Coulomb blockade prevents this. At first 
sight one might not expect that this repulsion could dra- 
matically modify the resonant current's dependence on 
the magnetic field. We find, however, that the resonant 
current is completely blocked for any value of the mag- 
netic flux except for integer multiples of the flux quantum 
$o = h/e. This striking effect has so far been overlooked 
in the literature, despite the current-switch in tunnel- 
coupled quantum dots, induced by coherent radiation of 
two microwaves, was proposed earlier 6]. 

In this Letter we show how such a switching effect can 
be seen in the exact analytical solution for certain sym- 
metric choices of the quantum dots' parameters. This 
solution can be found in the infinite bias limit, which 
we justify by comparing the finite- and infinite-bias cases 
for non-interacting electrons. We then present an ex- 
planation of the switching effect as a self-trapping phe- 
nomenon in non-linear systems. This explanation allows 
us to formulate general conditions for this phenomenon 
and clearly displays its generic nature. 

Consider a double dot (DD) connected in parallel to 
two reservoirs, as shown in Fig. [TJ For simplicity we 
consider spinlcss electrons. We also assume that each of 
the dots contains only one level, E\ and E 2 respectively. 
In the presence of a magnetic field, the system can be 



FIG. 1: Resonant tunneling through two parallel dots pierced 
by the magnetic flux $. 



described by the following tunneling Hamiltonian, 

H = H + H T + £/*4 d M + Ud Ui4d 2 . (1) 

H=l,2 

Here H = J2ki E kLa\ L a kL + E kR a\ R a kR \ and H T de- 
scribes the reservoirs and their coupling to the dots, 



H T = ^ \ + t^ R a\ R dA + H.c 



Lt.k 



(2) 



where fj, = 1, 2 and a\. L ^ R ^ are the creation operators for 

the electrons in the reservoirs while d\ 2 are the creation 
operators for the DD. The last term in Eq. {T]) describes 
the interdot repulsion. We assume that there is no direct 
transmission between the dots and that the couplings of 
the dots to the leads, t^^R), are independent of energy. 
In the absence of a magnetic field all couplings are real. 
In the presence of a magnetic flux however, the tunnel- 
ing amplitudes between the dots and the reservoirs are 
in general complex. We write t^um — l^LiR)^^^ , 
where is the coupling without the magnetic field. 

The phases around the closed circle are constrained to 
satisfy <j> 1L + tj>\ R - cj> 2 L - <t>2R = <j>, where <f> = 27r$/$ - 
Let the initial state of the system correspond to fill- 
ing the left and right reservoirs at zero temperature with 
electrons up to the Fermi energies /il and /jr, respec- 
tively (Fig. [T|). Consider first noninteracting electrons, 
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U = 0. In this case the problem can be solved exactly for 
any values of the bias voltage, pl — PR, and of the cou- 
plings to the leads Q • Indeed, the total wave function for 
the noninteracting electrons, \^(t)) — exp(— iHt)\^ (Q)) , 
can be written at all times as a product of single-electron 
wave functions, \^(t)) = T\ k , \ip k '(t)). Here \tpk'(t)) de- 
scribes a single electron initially occupying one of the 
levels Ek>L < PL, or Ek'R < pr in the left or right lead. 
It can be written as 

hM*)>=[ E b% k (t)al a + W*)<]|0), (3) 



k,a=L,R 



where b k j~ k (t) and b k 'fj,{t) are the probability amplitudes 
for finding the electron at the level EkL(R) m the left 
(right) reservoir, or at the level E^ inside the DD system. 
The total probability of finding the electron in the right 
lead is therefore P k ,(t) = J2k l & fc'fe(<)| 2 

Consider p^ > p R . Then the total average charge 
Q(t) accumulated in the right lead by time t is a sum 
of Pk> (t) over all electrons with initial energy within the 
potential bias, pr < Ek'L < Hl- The average current is 
I(t) = Q(t). (We adopt units where the electron charge 
e = 1). Let us take the continuum limit ^2 k — > J pdE k , 
where p is the density of state in the corresponding lead. 
Then the current can be written as 



I(t) = [ L p L I k ,(t)dE k , L , 

J PR 



(4) 



where Ik'(t) = dt \b k , k (t)\ 2 dE k R is a single electron 
current. 

Substituting Eq. ([3]) into the Schrodinger equation 
i\tp k '(t)) = H\ip k /(t)) we obtain the following equations 
for the amplitudes b(t), 

tt>k'k = E kL bi k + J2 (5a) 

11=1,2 

ibk'v = + £ (t^ L b L k , k + t; R b* k ) (5b) 

k 

ib k'k = E kRbk-k + E • ( 5c ) 

Replacing sum over the lead states by integrals, these 
equations can be solved analytically. Taking for simplic- 
ity tiL(R) = t2L(R) = ^l(r) we obtain for a single electron 
stationary current (t — > oo) the following result: 

h' = ^ R \t L \ 2 [ |M 2 + |/ 2 | 2 + 2Re(e- i */ 1 / 2 )] /D 2 (6) 

where T L(R) = 2<Kp L {R)\tL{R)\ 2 and 



fi.2 = (E k 'L-E 2A )T^(r-r^ 

D = (2E k , L -E 1 -E 2 + iT) 2 -e 2 + |r>| 2 . 
Here e = E 1 - E 2 , T = T L + T R , and T <t> = T L + T R e^. 



In the case of large bias, \pl,r — ^1,2 1 3> T, the inte- 
gration limits in Eq. ([4]) can be extended to infinity. As 
a result, we find for the total current I = I(<f>), 

e 2 +T L T R sm 2 

J (w = / o rx> ( 7 ) 

£ 2 + 4r L r fi sin 2 | 

where Iq = 2TlT r /T is the resonant current for non- 
interacting electrons in the absence of the magnetic 
filed. The (^-dependence in Eq.([7]) is an example of 
the Aharonov-Bohm effect; we illustrate this for finite 
and infinite bias in Fig. We find the infinite bias 
limit Eq. (|7|) is a very good approximation for a finite 
bias whenever the level is well inside the bias window, 

\pl,r - £1,2 1 > r. 

The entire treatment can in fact be simplified in the 
large bias limit by transforming Eqs. §5§ into Bloch-type 
master equations for the reduced density matrix of the 
DD system, a jy (t). Here j(j') = {0, 1, 2, 3} label the DD 
states in order: the empty DD, the first dot occupied, 
the second dot occupied, and both dots occupied. This 
density matrix is related to the amplitudes b k i ^(t) via Q 



bk'^t)bl,.,{t)pLdE k , L (8) 



for fi = \,2 and 000 = 1 - du - 022 — 033. Then multi- 
plying Eqs. ((5|) by b k ,At) and integrating over E k 'L gives 
the master equations 

(Too = -2r L cr o + T R (an + 022 + 012 + 021) (9a) 
oil = r L cr 00 - IVh + T R a 33 + ST*ai2 + STa 2 i (9b) 
022 = Tlctoo - ro-22 + r_R(T 3 3 + 5T*ai2 + STa2i (9c) 



CT33 = r L ((Tll + (T 2 2 - e l0 (Tl2 - 

CT12 = e^TLUoQ + ST (an + 022) 



J4> 



021) ~ 2T R a 33 (9d) 
TR.a 3 3 
- (ie + T)a u (9e) 

where ai 2 (t) = exp[i(tf>iR - </>2.r)]0"i2(*) and ST = 
(e^TL - T R )/2. The total current is given by U 

I(t) = T R [au(t) + 0-22 (t) + 2(733 (t) + 2Rea 12 (i)]. (10) 

Solving Eqs. I© in the steady-state limit (a —* for 
t — > 00) we easily reproduce Eq. ([7]). 

In contrast with the single-electron approach, the mas- 
ter equations ([9]) can be applied only for large bias. On 
the other hand, the master equations can accommodate 
the Coulomb blockade case (-^1,2 + U ~S> Pl) in the most 
simple way. Indeed, the Coulomb blockade merely leads 
one to exclude the states corresponding to a simultane- 
ous occupation of the two dots from Eqs. ([9]). As a result 
the master equations become 

ooo = -2T L a 0Q + T R (an + a 22 + a i2 + a 2 i) (11a) 

011 = r L o- o - r^o-ii - T R (ai2 + CT 2 l)/2 (lib) 
o"22 = r L cr o - T R a 2 2 - T R (ai2 + a 2 i)/2 (11c) 

012 = e 10 r L o- O o ^"( (Tl1 + (T22 ) _ + r «) 5 'i2 ( lld ) 
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where the total current is given by Eq. (| 10[) with 0-33 = 0. 
Equations (jTTJ) can in fact be rigorously derived by em- 
ploying the nonequilibrium Greens function in the weak 
coupling limit 0, 0, or directly from the many-body 
Schrodinger equation in the large bias limit [l(| by as- 
suming that max(T,T)/\nL,R. — -E1.2I <C 1, where T is 



11, hj 



Under this 



the temperature of the reservoirs 
assumption they are valid to all orders in the tunneling 
couplings. (In fact, the master equations approach has 
been already shown to be very useful for description of 
electron transport in coupled dots @,[l3, 14 1). 

Solving Eqs. (fTTj) in the steady-state limit we obtain 
for the total current 



m = io 



e 2 + I c ( 2T R sin 2 t - e sin < 



(12) 



where Ic = 2Ti,r R /(2Ti J + ^r) is the total current (with 
the Coulomb blockade) in the absence of the magnetic 
field. Comparing Eq. (JT2J) with Eq. © for e ^ we find 
that both currents display the Aharonov-Bohm oscilla- 
tions, Fig. [2] Nevertheless, Eq. (fT2|) shows an asymmetric 
behavior with respect to the magnetic flux, </> (the dashed 
curve in Fig. . It looks contradictory to the Onsager 
relation that locks the current peaks at 4> — 2im in any 
two-terminal linear transport [15| . Yet, this relation is 
not applicable in our case, corresponding to interacting 
system under finite bias voltage [4j . 

More strikingly is the current behavior in the inter- 
acting case for e — > 0. While the resonant current for 
the noninteracting electrons keeps oscillating with the 
magnetic field, Fig. [5^, it becomes non-analytic in <fi in 
the case of Coulomb blockade. Indeed, one finds from 
Eq. (fT2"j) that I — Ic for <f> = 2irn, where n — <5>/<l>o is an 
integer, but / = for any other value of $ (Fig. [5b). 

Such an unexpected "switching" behavior of the elec- 
tron current in the magnetic field can be explained in 
the following way. Let us disentangle the Coulomb block- 
ade and quantum interference effect, which interplay in a 
non-trivial way in the electron current. This can be done 
by defining new basis states of the DD, rf+ 10) -> 4l°>> 
chosen so that they do not interfere in the electron cur- 
rent. For instance, if the state a\\0) is not coupled to 
the right reservoir, i.e., t 2R — > t 2R = 0, then the cur- 
rent would flow only through the state d\\Q). Obviously, 
no interference between these two states would appear in 
the total current. 

Such a basis can be found for e = E\ — E 2 = when 
the DD Hamiltonian, Exd\di + E 2 d\d 2 , is invariant under 
SU(2) transformations. Then the unitary transformation 




Z2R 

* J.* 

2R L 1R 



di 
d 2 



(13) 




with N = (P 



1R 



i 2 .^) 1 / 2 , results in t 2 R = 0. 



FIG. 2: (a) Stationary current versus magnetic flux for the 
noninteracting DD with aligned and misaligned levels, where 
Fl = T_f? = r/2. Open circles show the result of numerical 
evaluation of Eqs. @, © for a finite voltage /j,l,r = Ei,2 ± 
5T; the curves show Eq.([7]), corresponding to infinite bias, 
(b) Stationary current for the same system in the Coulomb 
blockade regime, Eq (|12|l . 



Consider now the coupling of the state d) 2 \0) to the left 
lead. One obtains from Eq. ((2]) for t 2 L = t 2 L (</>) 

t 2L {<t>) = -e l ^-^\t lL i 2R e^ - t 2L i 1R )/M. (14) 

It follows from this expression that t 2 L = for (f> = 2mr 
provided that iii/^L = t\ R Jt 2 ^ or for 4> = {2n + l)ir if 
t\hlt 2 h = —i\R/i 2 R [H|]. Then the state d\\Q) decouples 
from both leads. This would not affect the resonant cur- 
rent for noninteracting electrons (U = 0), since the state 
d||0) is already decoupled from the right lead. 

In the case of Coulomb blockade, however, the coupling 
to the left lead becomes the point of crucial importance. 
Indeed, any discrete state coupled to an infinite reser- 
voir is going to be totally occupied, no matter how weak 
the coupling 0. Then the state d\\0), carrying the cur- 
rent, will be blocked by the Coulomb interdot repulsion. 
As a result, the total current vanishes. However, if the 
state d\\0) is decoupled from both leads, it remains unoc- 
cupied, so that the current can flow through the state 
d\\0). As shown above, this takes place precisely for 
t\hlt 2 h = iiin/t 2 R. If this condition is not fulfilled, 
the current is always zero, even for (f> = 2nn. 

It is clear from our explanation that the switching 
phenomenon disappears for e / 0, as can be seen 
from Eq. (|12| . Indeed, in this case the DD Hamil- 
tonian E^d/^dp is not invariant under SU(2) trans- 
formations. Therefore any unitary transformation that 
eliminates one link between the DD and the reservoirs, 
like Eq. ([±3]) . would generate direct coupling {t\ 2 = 
et 1R t 2R /[\t 1R \ 2 + \t 2R \ 2 }) between the states c? li2 |0) [H. 
As a result, the current is not interrupted by the Coulomb 
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FIG. 3: Transient behavior of the electric current for different 
values of the magnetic flux, obtained from Eqs. 



blockade, so that its behavior for e ^ would be similar 
to that in the noninteracting case, Fig. [2l 

The switching effect of the magnetic field on the elec- 
tric current takes place only in the stationary regime, 
t — > oo, where one of the states d| 2 |0) is fully occupied. 
However, it takes a certain transition time, ft r , during 
which the total current is not zero. One estimates from 
Eq. (JUJ) that tt t ~ l/r 2L , where T 2L = 2np L \t 2L \ 2 . Us- 
ing Eq. (HU) we find that t tr ~ r^ 1 ^/*) 2 for $ < <J> . 
Hence ttr becomes much longer than the usual relaxation 
time r _1 for a very small magnetic field (or in general 
when $ — » n$o). This is illustrated in Fig. 02 which 
shows the transient current I(t) as a function of t for dif- 
ferent values of <f>. One finds from this figure that the cur- 
rent always increases for small t. However, it eventually 
disappears for t 3> ttr- This explains the non-analyticity 
of I(4>) at <p = 27m. Indeed, at any finite time t there is 
no discontinuity at <j) = 2irn. It appears only in the limit 
of t — > oo because ttr — * oo for $ — > 0. 

As we demonstrated above, the switching effect be- 
comes very transparent in a particular basis of the DD 
states. Still, it is very surprising how such a basis emerges 
dynamically? Indeed, an electron from the left lead can 
enter the DD system in any of SU(2) equivalent super- 
positions of its states. Therefore there exists a probabil- 
ity for each electron to enter the DD in the superposi- 
tion that eliminates one of the links with the right lead. 
When it happens, the electron would be trapped in this 
state. Even if the probability of this event for one elec- 
tron is very small, the total number of electrons passing 
the DD goes to infinity for t — > oo. Therefore the trap- 
ping event is always realized for large enough time. In 
the presence of Coulomb blockade this would lead to the 
switching effect, as explained above. (A similar effect of 
the Coulomb blockade, leading to divergency in the shot- 
noise power has been discussed in different publications 
3 13, 12, 18 1 • However, this phenomenon is not related 
to vanishing of current, but rather to electron bunching 
leading to a system's bistability 17, 



Our interpretation of the switching effect allows us to 
determine the necessary conditions for its realization in 
real experiment. First, we need the total occupation of 
one of the states d\ 2 |0), decoupled with the right lead. 
This takes place only if the energy level is far from the 
corresponding Fermi energy, \hl — E\\ ^> Y. For instance, 
if E\ = [il, the occupation probability reaches only 1/2 
at t — ► oo. As the bias increases, however, other levels of 
the dots can enter into the bias window. In this case the 
transport would proceed through several levels. If the 
dots are identical, the same "rotation" (TIB"]) , applied to 
each pair of levels with the same energy, would result in a 
simultaneous decoupling of the corresponding "rotated" 
states from the right lead. Then the switching effect is 
expected to take place in this case as well, even if the 
inter-level spacing is very small. 

In real system, the surrounding environment will cause 
dephasing between the two dots due to fluctuations of the 
dots parameters, like energy levels, tunneling couplings, 
etc. Since a particular origin of these fluctuation is irrel- 
evant for evaluation of the corresponding dephasing rate 
[13] , one can model the environment by isolated fluctua- 
tors interacting with the system and vibrating its energy 
levels [Hf. Using such a model for the fluctuating envi- 
ronment one finds for the stationary current (for e = 
and T L =T R = T/2) Q, 



I 



l+/ c T d sin 2 (<^/2) ' 



(15) 



where I c = T/3 and the dephasing time. Therefore in 
order to observe the switching effect one requires I c Td 3> 
1/ sin 2 (0/2). For small 4> it is equivalent to 3> ttr- 

In general, for a weak coupling to the environment, 
the dephasing rate jd = 1/Td can be evaluated as jd ~ 
(Se) 2 S(Q), where Se is an average fluctuation of dots levels 
and S(u) is the corresponding spectral density [13, [20|. 
In the case of the thermal environment S(uj) oc T, where 
T is temperature [23]. Therefore by decreasing the tem- 
perature one can make the corresponding decoherence 
rate arbitrary small in order to reach the switching effect 
even for small values of the magnetic flux. 

As we explained above the experimental realization of 
the switching effect would require fulfillment of differ- 
ent conditions which should be met. The most essential 
of them are large bias and inter-dot Coulomb repulsion, 
which should exceed the bias. 

With respect to the bias, it would be hard to realize 
the large bias voltage (on scale of T) with one level inside 
the bias window. However, in the case of two identical 
dots, the switching effect is expected even if many levels 
are inside the bias. Therefore the large bias condition 
should not create essential experimental problem. 

A realization of the large inter-dot repulsion condition 
would represent a more complicated problem. Indeed, 
it implies that two dots are very close. This condition 
is still not met in present experiments [l[ . However, 
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it can be met by decreasing the dots size. An another 
way to achieve proximity of two parallel dots is to use 
different materials. For instance the quantum dots in 
graphene system can be a very promising set-up for an 
investigation of the switching effect (2f| . 
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